Spin splitting of X-related donor impurity states in an AlAs barrier. 
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We use magnetotunneling spectroscopy to observe the spin splitting of the ground state of an 
X-valley-related Si-donor impurity in an AlAs barrier. We determine the absolute magnitude of 
the effective Zeeman spin splitting factors of the impurity ground state to be g/= 2.2 ± 0.1. We 
also investigate the spatial form of the electron wave function of the donor ground state, which is 
anisotropic in the growth plane. 
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1. Introduction 



The effect of spin on electronic transport has recently 
attracted much interest both from an applied and from 
a fundamental point of view 0, 0- Resonant electron 
tunneling through the discrete states of self-assembled 
semiconductor quantum dots (QDs) and the physically 
similar bound states of impurities is a promising means 
of studying spin-resolved transport. Such tunneling ex- 
periments have been used to observe directly the spin 
splitting of the bound electronic states of shallow impu- 
rities within a GaAs quantum well 0, 0, or of elec- 
trons in InAs QDs H 0, @ , and to measure directly the 
absolute value of the ^-factor of these zero-dimensional 
states. In addition, they have been used to probe at a 
mesoscopic level the spin dependence of the local den- 
sity of states 0. Resonant tunneling studies have also 
been performed on single barrier GaAs/ AlAs/ GaAs het- 
erostructures. AlAs is an indirect gap material with the 
minima of the conduction band at the X points of the 
Brillouin zone, whereas in GaAs the minimum is at the 
r point. There have been several earlier studies of tun- 
neling through X- valley states [3 El El El 13 13 El , 
including investigations of the tunneling through donor 
states associated with the X conduction band minima 

[3 El El 13 M MM 

In this paper, we report the observation of Zeeman 
spin splitting of the ground state of an Si-donor impurity 
embedded in an AlAs tunnel barrier for the orientation 
in which the magnetic field is applied in the plane of the 
barrier, i.e. perpendicular to the direction of the elec- 
tron tunnel current. This state is associated with the 
anisotropic X-conduction band minima of AlAs. The Si 
atoms, which are located substitutionally on Al sites, dif- 
fuse during growth into the AlAs barrier from adjacent 
GaAs layers which are <5-doped with Si. We measure the 
effective g-factor of the zero-dimensional state and obtain 
the absolute values of the effective spin splitting factor 
components g of between 2.1 and 2.2. In addition, mag- 
netotunneling spectroscopy provides us with information 
about the spatial form of the wave function of an electron 
bound in the X- related donor state. Our measurements 



indicate that the wave function has a biaxial symmetry 
in the growth plane, with axes corresponding to the main 
crystallographic directions of the (001)- epilayers. 

Let us first briefly review the previously-studied prob- 
lem of tunneling through isolated donor impurities in the 
GaAs quantum well (QW) of large area GaAs/(AlGa)As 
double barrier resonant tunneling diodes. In this case 
the donors form localized (~10 nm) hydrogenic bound 
states associated with the T-conduction band minimum 
of the GaAs QW [24| . These states are located at an en- 
ergy of ~10 meV below the bottom of the lowest energy 
subband of the QW. Under an applied bias the tunnel 
current exhibits a rapid increase when the impurity state 
aligns with the Fermi level in the negatively-biased elec- 
tron emitter layer. In general, there are many impurities 
giving rise to multiple, overlapping steps in the current- 
voltage characteristics, and these multiple peaks can be 
resolved in the current-voltage characteristics of small 
area mesa samples 22]. Similarly, in our previous work, 
the resonant tunneling of electrons through individual X- 
related donor impurity states of a single, relatively thin, 
5 nm AlAs barrier (with an X conduction band quantum 
well) appeared as partially resolved fine structure in a 
broad resonance, associated with the ensemble of donors 
|2l|. This fine structure arises because the donors are 
located in different atomic planes of the AlAs and the 
spectrum of donor states is determined predominantly 
by the dependence of the binding energy on the position 
of the donor in the barrier. The influence of the random 
variations of the electrostatic potential on the energies of 
the donor impurity in this case is insignificant. 

In contrast, for the experiments described here, the 
donors are randomly located in a relatively thick, 11.2 
nm, barrier, so the influence of the random electrostatic 
potential is considerable. The essential role of the ran- 
dom variations of the electrostatic potential in this case 
is associated with the presence of the <5-doped layers near 
the barrier [25| and the slow dependence of the binding 
energy of donors on their position in the thick barrier |26| . 
As a result, resonant tunneling of electrons through the 
donor states gives rise to a series of sharp, well-resolved 
peaks in the /(V) curves. We ascribe each peak to tun- 
neling through a single or very small number of individual 
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donor states. This allows us to observe spin splitting of 
the donor resonances and to determine the g-factor of the 
zero-dimensional states directly. 

2. Samples 

A schematic diagram of our device is shown in Fig. 

I. The active part of our samples comprises a single 

II. 2 nm thick AlAs barrier which is sandwiched between 
two accumulation layers formed by two <5-doped layers 
with Si-concentration of 3 x 10 11 cm~ 2 , located at a 
distance of 2.8 nm from each side of the barrier. A two- 
dimensional electron gas (2DEG) forms in each accumu- 
lation layer at zero bias. The AlAs layer was not in- 
tentionally doped, but donor impurities are present in 
the AlAs due to diffusion of Si into the barrier from the 
(5-doped layers. The calculated T and X band profiles 
of the active part of our device at zero bias are shown 
in Fig.l. The heterostructures was grown by molecu- 
lar beam epitaxy on a (OOl)-oriented, Si doped n + -type 
GaAs wafer ( = 2x10 18 cm _3 ) at a tempera- 
ture of 550 C. The detailed layer composition of the 
heterostructure, in order of growth on the substrate, is 
as follows: a Si-doped, 0.5 /im-thick GaAs buffer layer 
(N d = 2xl0 18 cm~ 3 ), a 60 nm-thick GaAs layer (N d = 
3xl0 17 cm -3 ), a 21.6 nm-thick undoped GaAs layer; a 
5.6 nm-thick undoped Gao.9Alo.1As layer, a 28 nm-thick 
undoped GaAs layer, a Si (5-doped layer with concentra- 
tion of 3 x 10 11 cm -2 , a 2.8 nm-thick undoped GaAs 
layer, a 11.2 nm-thick AlAs barrier layer; a 2.8 nm-thick 
undoped GaAs layer, a Si J-doped layer with concen- 
tration of 3 x 10 11 cm~ 2 , a 28 nm-thick undoped GaAs 
layer, 5.6 nm-thick undoped Gao.9Alo.1As layer, a 21.6 
nm-thick undoped GaAs layer, a 60 nm-thick GaAs layer 
(N d = 3xl0 17 cm -3 ), and a 0.5 /im-thick, GaAs cap- 
layer (N d =2x 10 18 cm -3 ). Ohmic contacts were made by 
deposition and annealing of AuGe/Ni/Au layers. Mesa 
structures, with diameter between 50 /xm and 200 /im, 
were fabricated by chemical etching. 



3. Experimental 

Tunnel current measurements at constant applied volt- 
age with magnetic field B applied parallel to the current 
(i.e. perpendicular to the 2DEG) reveal Shubnikov-de- 
Haas (SdH)-like oscillations [27|. Close to zero applied 
bias, analysis of the SdH-like oscillations gives a value of 
n s = 3.27 x 10 11 cm -2 for the sheet density of the two 
2DEG layers. Figure 2(b) shows the low temperature 
(4.2K) current-voltage, I(V) , characteristics at low bias 
voltages for a typical device i which exhibit sharp peaks 
in the current over voltage range from 10 to 60 mV. This 
peak structure is observed to be sample-specific, but for 
a given sample it is exactly reproducible from one volt- 
age sweep to another. The peaks are reproducible even 
after thermal cycling of the sample, except for a small 



voltage shift (~ a few mV). We ascribe the peaks in cur- 
rent to single electron tunneling through individual, zero 
dimensional Si-donor states in the AlAs barrier. Similar 
features have been observed and reported previously in 
large area double barrier RTDs [3J, |5j, |28| and attributed 
to tunneling through individual tunneling channels due 
to zero- dimensional states. Increasing the voltage across 
the device moves the energy of the donor state relatively 
to the Fermi level of the 2DEG which acts as an emit- 
ter for the tunnelling electrons (Fig. 2(a)). Tunneling 
occurs as the donor state crosses the Fermi level of the 
2DEG, and stops when the donor state is brought below 
the 2DEG subband edge. 

Figure 2 show I(V) at 0T and 8 T with magnetic field 
oriented perpendicular to the current direction. In a 
magnetic field the ground state of an Si-donor impurity 
splits into two spin energy levels given by 

Esi = gi^BBm s (m s = ±1/2), (1) 

where gi is the g-factor of the Si-donor impurity. This 
opens up two separate channels for electrons from the 
2DEG to tunnel into, and we therefore see separate peaks 
in I(V) due to electrons tunneling through each of these 
spin energy levels. In a magnetic field applied perpendic- 
ular to the current (i.e. parallel to the 2DEG) the 2DEG 
emitter becomes partially spin polarized, due to energy 
splitting of the Fermi energies of the two spin species, as 
shown schematically in Figure 2(a). Due to the slow tun- 
neling rate from the 2DEG, the two spin species in the 
2DEG should be in thermal equilibrium, and so the chem- 
ical potential of each is the same. Therefore, there is an 
energy difference between the subband edge of the spin 
species, equal to the spin-splitting g2Df^BB, where g?D is 
the g-factor of electrons in the 2DEG emitter. Resonant 
tunneling occurs when an impurity spin level crosses the 
Fermi level of the 2DEG. We assume that spin is con- 
served during the tunnelling process. For each spin we 
see a separate onset of tunneling and the voltage dif- 
ference between the position of the onsets (AV pea k) is 
proportional to the energy difference, AE$i = gi^sB, 
obtained from Eq.(l). 

Figure 3(a) shows in detail the behavior of the first 
current peak of device ii at 4.2 K in different magnetic 
fields applied perpendicular to the current direction from 
B=0 T to 8 T, in steps of 0.5 T. Figure 3(b) shows the 
voltage separation between the corresponding two spin- 
split peaks which increases linearly with magnetic field 
strength, as expected for a Zeeman effect. Because of 
the finite widths of the current peaks, it is not possible 
to resolve the splitting for magnetic fields less than 5 T. 
The best fit line to the data closely intersects AV=0 at 
B=0 T and has a slope gifis / /, where [Ib is the Bohr 
magneton, gi is the effective gyromagnetic ratio of the 
impurity with the magnetic field perpendicular to the 
current (i.e. perpendicular to the growth direction of the 
quantum well), and / is the so-called electrostatic lever- 
age factor. The temperature dependence of the current 
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onset allows us to determine the electrostatic leverage 
factor / 0, EJ . Figure 4 shows the temperature depen- 
dence of the onset of a typical resonant current features. 
We deduce a value of 0.44 for the electrostatic leverage 
factor by fitting the Fermi-Dirac function to the form of 
the measured low bias onset of the peak in current at 
various temperatures, using the procedure described in 
reference @. A similar value of / is obtained from self- 
consistent Poisson-Schrodinger calculations: these indi- 
cate that, over the bias range of interest (0-100 mV), the 
leverage factor / for an electron tunneling from the emit- 
ter into an impurity located at the center of AlAs barrier 
varies slightly from 0.44 to 0.42 eV per volt of applied 
bias. The uncertainty in our value of g is determined by 
the error in the leverage factor. Note also from Figure 
4 that the resonant peak in I(V) is strongly enhanced 
as the temperature is reduced. The enhancement, which 
may be related to a many-body Fermi energy singularity 
effect [23, tends to improve the resolution of the spin- 
split peaks. 

We have measured the splitting of several peaks in 
I(V) and find that the different impurity-related peaks 
give values of gi in the range from 2.1 to 2.22. Our value 
of the (/-factor of the X-related impurity states in AlAs 
is of a larger absolute value than reported in another 
tunneling experiment |20|. where g—OSA. However, for 
the experiment described here, the donors are located 
in a relatively thick, 11.2 nm, AlAs layer, whereas the 
localized state investigated in Ref. [2(j was embedded 
in a narrow 2 nm AlAs barrier. This value of (7=0.34 
is quite different from that for the X-valley electrons in 
bulk AlAs. The (/-factor for electrons in bulk AlAs ex- 
pected from theoretical calculations is 1.9 [3(J, and the 
g-factor of electrons in bulk Al0.sGa0.2As has been mea- 
sured by electron paramagnetic resonance to be 1.96 |3lj . 
Also, van Kesteren et al. have reported a value of « 1.97 
for electrons in AlAs QWs based on optically detected 
magnetic resonance experiments on AlAs-GaAs super- 
lattices [s^. This difference in values may be due to the 
complex nature of the X-related donor impurity states in 
AlAs barrier J2jj. In low-dimensional heterostructures, 
it is known that the value of the g factor can be modified 
from its bulk value, due to quantum confinement effects 
and because the electron wave function contains contri- 
butions from the different materials which make up the 
structure Calculations show that g is a strong 

function of the quantum well width |35j . and the rather 
wide X minima quantum well of our AlAs barrier gives a 
g-factor of the X-donor in our experiment which is close 
to the g- value of bulk AlAs g s=s2, since the modification 
of the band structure due to the quantum confinement is 
fairly small. 

We also studied the g-factor dependence for magnetic 
field applied along different crystallographic axes in the 
(001) growth plane. In contrast to the case of single- 
electron tunneling through the localized states of InAs 
quantum dots [8j , to within the resolution of our experi- 
ment, the spin splitting of the X-related donor impurity 



in AlAs was isotropic with respect to the angle of the 
in-plane magnetic field. 

We now consider the magnetic field dependence of the 
amplitude of the tunnel current through the X-related 
donor impurities as a function of magnetic field B ap- 
plied perpendicular to the direction of tunneling. We 
attribute the general fall in amplitude of both spin-split 
components with increasing B (see Fig. 5(a)) to a well- 
established effect that can be understood in term of a 
single-particle model for electron tunneling in the pres- 
ence of magnetic field |U 0, . Let a, (3 and Z indi- 
cate, respectively, the direction of B, the direction nor- 
mal to B in the growth plane (X, Y), and the normal to 
the tunnel barrier, respectively. When an electron tun- 
nels from the emitter accumulation layer into the impu- 
rity state in the barrier, it acquires an additional in-plane 
momentum given by 

kp = eBAs/h, (2) 

where As is the effective distance tunnelled along Z{ 
~ 8 nm for our device ). This gives the increased momen- 
tum along 0, which is acquired by the tunnelling electron 
due to the action of the Lorentz force. 

The applied voltage allows us to tune resonantly to the 
energy of a particular impurity state. Thus, by measur- 
ing the variation of the tunnel current with B, we can 
determine the size of the matrix element that governs 
the quantum transition of an electron as it tunnels from 
a state in the emitter layer into an impurity state. 

In order to analyse the results of our experiment, we 
express the tunneling matrix element M in terms of 
the Fourier transforms 3>j(/) (k) of the conventional real 
space wave functions, according to the relation M = 
J. <&i{k — kp)<f>f(k)dk, and express the tunneling current 

as I ~ \M\ 2 |23,l23- Here tne subscripts i and / in- 
dicate the initial (emitter) and final (Si-impurity) states 
of the tunnel transition. Relative to the strong spatial 
confinement in the zero dimensional impurity state, the 
initial state in the emitter is essentially unconfined - i.e. 
it behaves like a free particle. Hence, in /c-space 3>i(k) 
corresponds to a sharply peaked function with a finite 
value only close to k— 0. Since the tunnel current is 
given by the square of the matrix element involving <I>i (k) 
and the Fourier transform of the Si-donor wave function, 
&Si(k) , the narrow spread of k for <F;(fc) allows us to 
investigate the form of &si(k) by varying B and hence 
k, according to equation 2. Thus by plotting 1(B) for 
a particular direction of B we can measure the depen- 
dence of |<E>si(fc)| 2 along the k -direction perpendicular 
to B. Then, by rotating B in the plane (X, Y) and mak- 
ing a series of measurements of 1(B) with B set at regular 
intervals (A8 ~ 15 ° ) of the rotation angle 9 , we obtain 
a full spatial profile of \$>si(kx> ky)\ . This represents 
the projection in /c-space of the probability density of a 
given impurity electronic state peak. 

Typical experimental data for sample it, of the varia- 
tion of the current at peak A with the direction of mag- 
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netic field at 8 T, are shown as polar plots in Fig. 5 (b) and 
(c). A maximum current modulation AI / 1 of about ~ 
23% is observed, with clear, twofold anisotropy observed 
for the two split peaks A - and A + . The anisotropy of 
each observed peaks has a similar magnitude and orien- 
tation. In Fig. 5, 0° corresponds to the [110] direction, so 
the principal axes for the anisotropy are oriented along 
[100] and [010] directions. This result shows that the 
wave function shape of the X-related donor impurity in 
AlAs barrier is anisotropic in the growth plane, with the 
wave function probability density elongated along the di- 
rection [100] in real space. This is in contrast to the 
case of Si donor states in a GaAs quantum well, where 
the electron wave function has circular symmetry in the 
growth plane, as expected for Is donor ground state |37| . 
We suggest that this anisotropy may be related to the 
anisotropy of the effective mass of the electrons in the 
X conduction band minima of AlAs, though this point 



requires further theoretical analysis. 

In conclusion, we have used magnetotunneling spec- 
troscopy to observe the spin splitting of the ground state 
of Si donor impurities in an AlAs tunnel barrier. These 
states are associated with the X-conduction band min- 
ima of AlAs. We determine the absolute magnitude of the 
anisotropic effective magnetic spin splitting factors, g, for 
these states to be 2.1 ± 0.1. In addition, we use magneto- 
tunneling spectroscopy to investigate the spatial form of 
the wave function of the X-related donor impurity. The 
wave function of electrons bound to an X-related donor 
has a biaxial symmetry in the growth plane, with axes 
corresponding to the main crystallographic directions. 
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FIG. 1: Calculated conduction band profile of the active part 
of the tunnel structure at zero applied voltage. The figure 
shows the positions of the Fermi level Ef, the quantised 
GaAs accumulation layer subbands Eql and Eon, and the 
size-quantized levels of the Xz and Xxy subbands in the 
AlAs barrier. The energy positions of the X-impurity related 
levels in AlAs are also shown. 
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FIG. 2: (a) Schematic diagram of the spin splitting of the 
Si donor state in the AlAs and partial spin polarization of 
the 2DEG in a magnetic field applied perpendicular to the 
current. The effect of applying a voltage across the device 
is move the energy donor energy levels down relative to the 
Fermi level of the 2DEG. (b) /(V) characteristics at 4.2 K 
of sample i at T and 8 T for magnetic field applied per- 
pendicular to the current. The curves are offset for clarity. 
The characteristics show sharp peaks in the current due to 
tunneling through discrete X-impurity states, each of which 
split in a applied magnetic field. 
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FIG. 3: (a) Evolution of the spin-split peaks in I(V) the 
sample ii at 4.2 K in various magnetic fields perpendicular to 
the current from B=0 T to 8 T, in steps of 0.5 T. (b) The 
measured spin splitting versus magnetic field. The dashed 
lines are linear fits to the data. 




FIG. 4: I(V) characteristics of the first current peak of the 
typical samples at different temperatures showing the Fermi- 
level broadening. 
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FIG. 5: (a) Amplitude dependence of the peak A of sample h 
vs magnetic field B applied parallel to a [110] direction in the 
plane of the quantum well. Polar plot of the change in peak 
current vs in-plane magnetic field direction for the peaks A+ 
(b) and A- (c) at 8 T. 



